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ABSTRACT 

We investigate the serendipitous X-ray source population revealed in XMM-Newton ob- 
servations targeted in the Galactic Plane within the region 315°< I < 45° and |6| <2.5°. 
Our study focuses on a sample of 2204 X-ray sources at intermediate to faint fluxes, which 
were detected in a total of 1 16 XMM-Newton fields and are listed in the 2XMMi catalogue. 
We characterise each source as spectrally soft or hard on the basis of whether the bulk of the 
recorded counts have energies below or above 2 keV and find that the sample divides roughly 
equally (56%:44%) into these soft and hard categories. The X-ray spectral form underlying 
the soft sources may be represented as either a power-law continuum with F ^ 2.5 or a ther- 
mal spectrum with kT ^ 0.5 keV, with Nh ranging from 10^0^22 (.m"^. For the hard sources, 
a significantly harder continuum form is likely, i.e., F ^ 1, with Nh = 10^^^^^ cm^^. For 
^50% of the hard sources, the inferred column density is commensurate with the total Galac- 
tic hne-of-sight value; many of these sources will be located at significant distances across the 
Galaxy implying a hard band luminosity Lx > 10^^ erg s^^, whereas some will be extra- 
galactic interlopers. A high fraction (^90%) of the soft sources have potential near-infrared 
(2MASS and/or UKIDSS) counterparts inside their error circles, consistent with the domi- 
nant soft X-ray source population being relatively nearby coronally-active stars. These stellar 
counterparts are generally brighter than J=16, a brightness cutoff which corresponds to the 
saturation of the X-ray coronal emission at Lx — 10^^ L^oi- In contrast, the success rate in 
finding likely infrared counterparts to the hard X-ray sample is no more than « 15% down 
to J=16 and w 25% down to J=20, set against a rapidly rising chance coincidence rate. The 
make-up of the hard X-ray source population, in terms of the known classes of accreting and 
non-accreting systems, remains uncertain. 
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1 INTRODUCTION 

The brightest sources discovered in the first all-sky surveys con- 
d ucted at X-ray wavele ngths (e.g., Uhuru, Forman et al. 191^ Ariel 
K lWarwicketaI.1 19811) , were found to be X-ray luminous close bi- 
nary systems powered by the accretion of matter onto a compact 
object. These sources, with intrinsic X-ray luminosities (Lx) typ- 
ically in the range IQ'^''"'^* erg s~^ in the 2-10 keV band, can be 
classed either as low-mass or high-mass X-ray binaries (LMXBs or 
HMXBs) depending on the nature of the non-degenerate star, the 
companion objects being either neutron stars or black holes. X-ray 
catalogues of the time also contained a few examples of other types 
of Galactic X-ray source, including X-ray bright supernova rem- 
nants, cataclysmic variables (CVs) and coronally-active binaries. 
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such as RSCVn systems, albeit with inferred X-ray luminosities 
typically much less than 10'^® erg s~^. 

Later more sensitive X-ray surveys utilising imaging instru- 
ments operating in a somewhat softer spectral regime (e .g., Ein- 
stein, Hertz & Grindlay 1984; ROSAT. Woees et alJI 19991) showed 
that the X-ray sky at low Galactic latitude is crowded, with nearby 
coronally-active stars and binaries dom inating the source statis- 
tics in the soft « 2 keV) X-ray band jMotch et af]|l997h . Sub- 
sequen tly, an imagine surve y in the Galactic plane carried out by 
ASCA jSugizaki et al]|200ll) provided a first detailed view of the 
general population of faint Galactic X-ray sources in the 2-10 keV 
band. Since the impact of X-ray absorption by interstellar gas is 
greatly diminished above ^ 2 keV, the ASCA survey was able to de- 
tect hard-spectrum sources with X-ray luminosities markedly lower 
than those of the classical Galactic X-ray binaries across a signifi- 
cant fraction of the inner Galaxy. 

More recently, the improved sensitivity, spatial resolution and 
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energy range afforded by the Chandra and XMM-Newton ob- 
servatories has provided the opportunity to revisit the issue of 
the faint Galactic source populations in both the soft and hard 
X-ray bands. For both missions, the fields-of-view of the on- 
board X-ray cameras are such that in a typical low-latitude point- 
ing, in addition to the primary target, many tens of serendipi- 
tous sources are seen. Both missions have also invested observ- 
ing time in studying specific regions of the Galactic plane, most 
notably the Galactic centre, through both deep pen c il-beam ob- 
serva t ions I Ebisawa et alj 200ll: Ebisawa et al. 20051: Muno et al.l 



l2003l : lMuno et alj|2004l : iHong et al.ll2009l : iRevnivtsev et aLlllOOsI) 
and mini-surveys in which wider angle coverage is ach ieved by the 
mosa i cing of multiple (relatively short) observatio ns jWang et al.l 
2OO2I: iHands et al.ll2004l : Iwiinands etZll2006 : Muno et al.ll2006l : 



Muno et al.ll2009h . The outcome is that 10 years post-launch both 



the Chandra and XMM-Newton archives contain a wealth of data 
relevant to Galactic X-ray sources at intermediate to faint flux lev- 
els, spanning a wide range of intrinsic luminosity. 

Recent results from Chandra and XMM-Newton demonstrate 
the potential of Galactic X-ray surveys, at readily accessible sensi- 
tivity limits, to detect a wide variety of source types. For instance, 
coronally-a ctive binaries can b e detected at distances up to 1 kpc 
or beyond ( lHerentetal]|2006h and young stellar objects (YSOs) 
can be unveiled in regions of current star-formation. In the case of 
the latter, evolved protostars and TTauri stars with extreme coro- 
nal emission and hard X-ray spectra (kT~l-4 keV) can be detected 
in dense mo lecular clouds despite the la rge line-of-sight column 
density (e.g. jFeigelson & Montmerlj 1999 1 Isolated Neutron Stars 
(ISNs), such as those discovered in the ROSAT survey (Haberl 
2007), are radio-quiet objects located at distances of no more than 
a few hundred parsecs that display soft thermal X-ray spectra. CVs 
constitute a source class that may, potentially, be found in large 
numbers in sensitive Galactic X-ray surveys. In particular, it has 
been proposed that intermediate polars may account for a large 
fraction of th e low Lx sour ces which reside in the Galactic Cen- 
tre region dMuno et alj2006l) and through their integrated emission 
may account fo r a significant fraction of the hard Galactic X-ray 
Ridge emission jSazonov et alj|200d IRevnivtsev et alj|2006h . Rel- 
atively quiescent X-ray binaries with either low-mass (e.g., black 
hole transients) or high-mass {e.g., Be star) secondaries may also 
make a non-negligible contribution to the source statistics. 

However, in truth, our knowledge of the makeup of the Galac- 
tic X-ray source population at relatively faint levels is quite limited. 
This is certainly the case in the 2-10 keV band where, in principle, 
the visible volume extends to the edge of the Galaxy. In order to 
better define the various populations in terms of their space den- 
sity, scale height and luminosity function, we need to characterise 
and, where possible, identify much larger samples of sources than 
are currently available. More comprehensively characterised source 
samples might also reveal how X-ray sources map onto structures 
such as the Galactic spiral arms, the thin and thick disc, the Galactic 
bulge and the mass concentration within ~ 100 pc of the Galactic 
Centre. A range of astrophysical issues, for example, relating to the 
formation and evolution of close binary systems and accretion at 
low mass-transfer rates, might also be addressed. 

In the above context, the Chandra Multiwavelength Plane 
(ChaMPlane) Survey is aiming at a systematic analysis of low- 
latitude fields with the objective of measuring or constraining the 
populations of low-luminosity ( Lx ^ 10'^^ erg s^^) accreting 
white dwarfs, neutron stars, and stellar mass bla ck holes in the 
Galactic plane and bulge (see lGrindlav et al1l2005l for full details). 
The programme utilises Chandra X-ray data in combination with 



follow-up optical and IR photometric and spectroscopic observa- 
tions, so as to maximise the number of identified sources and ex- 
plore the populations thereby revealed. Recent publications from 
the ChaMPlane program me include a stud y of Chandra fields in 
the Galactic anti-centre taong et alj|2005l) and the Galactic cen- 
tre and bulg e ('Lavcock et al.'l2005l: lKoenig et af]|2008l : lHong et al.1 
I2OO9I : Ivan d en Berg et al. 2009). 

In the case of XMM-Ne wton, the XMM-N ewton Galactic 
Plane Survey (hereafter XGPS: iHands et alj|2004l) has sampled a 
flux range which bridges the gap between the rel atively shallow 
ASCA Galactic Plane Survey jSugizaki et al.ll200l]) and the sensi- 
tivity limits reached in deep Chandra poi ntings {e.g., Ebisawa et al.l 
l200lh . In a paper reporting the XGPS, iHands et al.l ( 12004 ') dis- 
cussed the results from a programme including 22 pointings which 
cover a region of approximately 3 square degrees between 19°- 
22° in Galactic longitude and ±0.6° in latitude. Subsequent opti- 
cal follow-up observations of a representative sample of the bright- 
est identified low-latitude hard band sources detected in the XGPS 
have recently been presented by Motch et alj ( 1201 0). 
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by carrying out a systematic investigation of 
the X-ray source population seen at intermediate to faint fluxes 
in XMM-Newton pointings targeted at the Galactic plane. More 
specifically, we consider observations encompassing a narrow strip 
of the Plane towards the central quadrant of the Galaxy. In the next 
section of this paper, we give details of the set of XMM-Newton ob- 
servations which provide the basis of the study. In Section 3, we de- 
scribe the selection criteria we employ to extract a clean "serendip- 
itous" X-r ay source sample f rom these fields, using the 2XMMi 
catalogue dWatson et al.l2009l) as the input database. Section 4 goes 
on to investigate various properties of the sample with a focus on 
the available X-ray spectral information. In section 5 we present 
an investigation of the likely near-infrared counterparts based on 
a cross-correlation of the X-ray sample with both the 2MASS and 
UKIDSS surveys. This leads on to a discussion of the nature of the 
soft- and hard-source populations in section 6, followed by a brief 
summary of our conclusions of this first paper (Paper I). In sub- 
sequent papers (in preparation), we will explore the average X-ray 
spectral properties of a subset of relatively bright sources drawn 
from our source sample (Paper II), the log N - log S curves for both 
the soft and hard source samples (Paper III), and the broad-band 
colours of likely counterparts to the X-ray sources (Paper IV). 



2 THE XMM-NEWTON GALACTIC PLANE DATABASE 

In this paper we utilise observations drawn from the XMM-Newton 
public data archive targeted at positions along the Galactic plane 
within the central quadrant of the Galaxy. More specifically we use 
observations with pointings in the region bounded by 315° < I < 
45° and |fo| <2.5°. Our preliminary list of observatio ns comprised 
those used in constructing the 2XMMi catalogue IWatson et al.l 
I2OO9) - see Fig. [T] However, given our focus on the serendipi- 
tous source content of the XMM-Newton fields, we excluded those 
fields dominated by a very bright source, which in most cases was 
the target source. We also excluded a number of observations oth- 
erwise dominated by the target, such as those containing nearby 
star clusters or bright extended emission arising from supernova 
remnants. In some instances we excised the region within the field 
of view dominated by the bright emission but retained the rest of 
the observation. Observations which significantly overlapped with 
other deeper observations were also removed, although fields with 
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fractional overlaps giving mosaiced coverage of particular regions 
were retained. A couple of observations with relatively short net 
exposures affording only poor sensitivity were also discarded. 

This filtering process resulted in a final list of 116 observations 
which comprise the present survey. Details of the set of selected 
observations are given in Table [T] and the corresponding distribu- 
tion of pointing directions along the Galactic plane is illustrated in 
Fig-El In brief, Table[T]provides the following information for each 
field: column 1 - the XMM-Newton observation ID; columns 2-3 - 
the nominal observation pointing position in Galactic coordinates 
(l,b); columns 4-6 - net exposure time in the three EPIC cameras; 
columns 7-9 - the filter deployed in each EPIC camera; column 10 
- the fraction of the total field of view retained in those fields in 
which a sub-region was excluded; columns 11-13 - the number of 
sources detected in the field which pass the selection criteria de- 
fined in §3 (both the total number and the split between the soft and 
hard assignments are given). 

Most of the observations wer e performed with the three EPIC 
came ras (pn, MOS-1 and MOS-2; [striider et alj|200ll ; [Turner etal] 
l200lh in operation in a Full-Frame Window mode for which all 
the pixels of all the CCDs are read out and the full field-of-view 
is covered. For five of the observations no suitable pn data were 
available. Although, in most cases, the Medium filter was selected, 
there were a few occasions when either the Thinl or Thick filters 
were deployed (see Table [TJ. The distribution of the effective ex- 
posure time, calculated as the average of the individual pn, MOS 1 
and M0S2 exposures, with a double weighting for the more sen- 
sitive pn instrument, is illustrated in Fig.|2l the minimum effective 
exposure was 1.5 ks and the maximum 71.3 ks. The total sky area 
covered by the set of survey fields was 22.5 sq.deg. 



3 DEFINING A SERENDIPITOUS SOURCE SAMPLE 

As the starting point for the compilation of a serendipitous X- 
ray source sample, we extracted a complete list of the sources 
detected in the 116 designated fields from the S econd XMM- 
Newt on Serendipitous Source Catalogue (2XMMi; I Watson et al.l 
l2009l) . The 2XMMi catalogue provides extensive details of the 
source properties including the net count rate and detection max- 
imum likelihood of each source in five separate energy bands: 
Band 1 (0.2-0.5 keV); Band 2 (0.5-1.0 keV); Band 3 (1.0-2.0 keV); 
Band 4 (2.0-4.5 keV): and Band 5 (4.5-12.0 keV). This information 
is provided for each of the three EPIC instruments. 

A preliminary step involved dealing with the relatively small 
number of observations in which instrument filters other than the 
medium filter were deployed. In these cases the recorded count 
rates measured in each camera and in each energy band were scaled 
to equivalent values for the medium filter, where the scale factors 
applied were derived from the r atio of the appropria te energy con- 
version factors (ECF) quoted in lMateos et al.l ( l200^ n. 

In our study of intermediate to faint Galactic Plane sources, 
it has proven convenient to compress the available spectral infor- 
mation into two bands, a soft band representing the combination 
of Bands 2 and 3, encompassing the energy range 0.5-2 keV, and 
a hard hand based on Bands 4 and 5, corresponding to 2-12 keV. 
The 2XMMi Band 1 information has not been used in this analy- 
sis. The count rates measured in each EPIC camera for the soft and 
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Figure 1. The spatial distribution of the XMM-Newton observations that 
fall within the Galactic plane region bounded by |(| < 45° and |6| <2.5°. 
Each colour dot co iTesponds to a fie ld used in the construction of the 
2XMMi catalogue jWatson et alj|200 9), with the colour coding red-to- 
green-to-blue roughly scaling as the logarithm of the accumulated exposure. 
The 1 16 observations included in the present study are circled in black. 
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Figure 2. The distribution of the effective exposure time across the set of 
116 observations, calculated as the average of the individual pn, MOS 1 and 
MOS2 exposures with a double weighting for the more sensitive pn instru- 
ment. 



^ These ECFs are calculated assuming a broad-band source spectrum, char- 
acterised as a power-law continuum with photon index F = 1.7 absorbed by 
a line-of-sight column density Nfj = 3x lO'^" cm~^. 



4 R.S. Warwick, D. Perez-Ramirez and K. Byckling 

Table 1. Details of the 116 XMM-JVewton observations which comprise the current Galactic Plane survey 



Observation 
ID 


Pointing 
I b 


Exposure (ks) 
pn MOSl MOS2 


pn 


Fiher 
MOSl 


M0S2 


FoV 
Fraction 


total 


Sources 
soft 


hard 


0400010901 


318.9 


+0.5 


9.5 


11.6 


11.6 


Med 


Med 


Med 




13 


12 


\ 


040S'^Q04-ni 


323.4 


-1.7 


12.4 


16.6 


16.3 


Med 


Med 


Med 


97 


13 




5 


OAOOi^QO 1 1 


325.9 


-1.8 


15.4 


21.6 


23.2 


Thinl 


Thinl 


Thinl 


88 


27 


19 


c 
o 


090^01 01 01 


327.2 


-0.1 


7.1 


8.7 


8.7 


iVlCU 


iVlCU 


Thinl 


Q5 


16 


10 




1 zl^'^R^ 1 1 


328.4 


+0.2 


0.0 


50.0 


50.2 




Thick 


Thick 





^4 


26 


c 
o 


VJH-UtHJJU 1 U 1 


331.7 


-0.6 


24.7 


31.6 


31.6 


Med 


Med 


Thinl 




f\\ 

\ji 


47 


14 


040/^7 09 01 


333.8 


+0.1 


6.0 


7.6 


7.6 


iVlCU 


iVlCU 


iVlCU 




1 R 


1 % 




040^9 R09 01 


334.7 


+0.0 


7.5 


9.9 


11.1 


iVlCU 


iVlCU 


M(»/l 


ft 89 


zo 




7 
1 


0901 000401 


335.6 


-0.4 


18.1 


21.0 


21.3 


Thinl 


Thinl 


Thinl 


ft Q8 

0.70 


^9 

jZ 


90 
zu 


12 


0099140101 


336.4 


-0.2 


4.3 


13.8 


14.7 


Med 


Med 


Med 


0.97 


g 


5 


3 


0904^00^^01 


337.1 


+0.0 


30.7 


32.6 


32.6 


Thinl 


Thinl 


Thinl 


ft 88 


9R 

Zo 


1 s 

ij 


1 


0^^07170901 


337.4 


-0.3 


60.0 


81.2 


84.0 


Med 


Med 


Med 




70 


48 


22 


0^0^100101 


337.8 


-0.1 


38.3 


45.3 


45.7 


Med 


Med 


Med 


0.97 


51 


37 


14 


0900Q001 1 


343.1 


-0.1 


23.9 


27.3 


27.6 


Med 


Med 


Med 




9R 

Zo 


1 s 


1 '\ 


04.0(^7^9^01 


343.9 


+0.0 


2.7 


3.5 


4.3 


Med 


Med 


Med 




9 


3 




01 19460901 


345.0 


+2.5 


20.9 


27.7 


27.9 


Med 


Med 


Med 




26 


12 


14 


04067'iO'^OI 


345.3 


+0.0 


6.1 


7.6 


7.6 


Med 


Med 


Med 




18 


13 


5 


0144080101 


345.7 


-0.2 


10.9 


16.9 


16.9 


Med 


Med 


Med 




28 


18 


10 


0401 Q601 01 

O^O 1 !7UO 1 1 


348.9 


-0.5 


15.5 


19.2 


20.3 


Med 


Med 


Med 




28 


18 


10 


niO'^961901 


350.5 


-0.5 


4.3 


6.7 


6.7 


Med 


Med 


Med 






5 


3 


01 19901401 


354.1 


+0.1 


4.6 


6.3 


6.4 


Med 


Med 


Med 


77 


7 


3 


4 


0^017^0101 


357.6 


-0.1 


59.8 


61.5 


61.5 


Thick 


Thick 


Thick 


0.95 


67 


40 


27 


010'^9^^1 ^01 


358.6 


-1.0 


5.0 


7.6 


7.6 


Med 


Med 


Med 




10 


2 


g 


090(SQQ0 1 1 


358.6 


-2.2 


7.9 


13.5 


13.5 


Med 


Med 


Med 




43 


28 


15 


040^1 S 80901 


358.8 


-0.5 


14.1 


18.4 


19.1 


Med 


Med 


Med 




21 


9 


12 


0400*^401 01 


359.1 


-0.5 


17.2 


20.2 


22.3 


Med 


Med 


Med 


0.95 


25 


10 


15 


n'^049901 01 


359.3 


+0.9 


0.0 


7.5 


7.7 




Med 


Med 


0.98 


5 


4 


\ 


01 1 9Q70701 


359.6 


+0.0 


19.4 


23.4 


23.4 


iVlCU 


iVlCU 


iVlCU 




jy 


1'^ 


Afx 

H-D 


0'^09RR4'^01 


359.7 


+0.3 


5.0 


6.5 


6.5 


iVlCU 


iVlCU 


iVlCU 




IQ 

i.y 


o 


1^ 


01 19Q71001 


359.9 


-0.3 


11.3 


11.3 


11.5 


Thick 


Med 


iVlCU 




zo 


1 < 


w 


1 1 9Q71 QO 1 


0.0 


+0.3 


4.2 


7.5 


7.2 


Med 


Med 


Med 




\ \ 




Q 
o 


0^071 imni 

UjU / i iUiUi 


0.0 


-1.6 


13.3 


21.2 


21.9 


ivieu 


iVicU 


iVlcCl 




OA 


1 4 


1 n 




0.3 


-2.2 


17.0 


20.5 


20.5 


iviea 


iviea 


iviea 




AO 


ZV 


1 1 
1 J 


01 1 9071 SOI 


0.6 


+0.0 


6.9 


9.3 


9.3 


ivieu 


ivieci 


ivieu 




14 


7 


7 
1 


01 1 9Q70901 


0.9 


+0.1 


12.9 


17.4 


17.4 


iVlCU 


iVlCU 


M(»/1 


ft Q8 


^R 

JO 


1 < 

1 J 


Zj 


0'^0S8'^0701 


0.9 


+0.8 


3.0 


4.0 


4.4 


Med 


Med 


Med 




3 


\ 


2 


090^9401 01 


1.0 


-0.1 


15.1 


32.8 


34.3 


iVlCU 


iVlCU 


iVlCU 




40 


ID 


24 


n'^09^70101 


1.0 


-1.9 


27.1 


31.1 


31.1 


Med 


Med 


Med 




60 


34 


26 


01 4R000S01 


3.3 


-0.3 


6.5 


8.8 


8.7 


Thinl 


Thinl 


Thinl 




iu 


Q 
o 


Q 
o 


0^01 SS1001 


6.1 


+ 1.2 


14.1 


15.9 


14.8 


Med 


Med 


Med 




17 


9 






6.1 


-1.9 


18.2 


21.3 


21.2 


Thin 1 

i nini 


Thin 1 


Thin 1 

1 nin 1 




JO 


1 
Zl 


1 s 


01 '^5749^01 1 


6.7 


+0.3 


3.5 


4.1 


4.2 


Med 


Med 


Med 


81 


4 


3 


I 


01 '^S749S01 1 


6.8 


+0.0 


2.3 


2.8 


2.9 


ivieu 


iVlcU 


ivieu 


A 71 


■2 
J 


z 


1 
i 


01 '^5749601 1 


7.0 


+0.3 


4.0 


4.9 


4.9 


Med 


Med 


Med 




10 


g 


2 


01 '^S749R01 1 

\JLJJ /tZoUl ' 


7.2 


+0.0 


4.7 


5.8 


6.0 


iVlCU 


iVlCU 


iVlCU 




17 


ID 


I 


01 ■^S749Q01 1 


7.3 


+0.3 


4.4 


6.1 


6.3 


Med 


Med 


Med 




12 


10 


2 


01'^S74'^10lt 

\JLJJ /tj iUl ' 


7.5 


+0.0 


4.9 


6.6 


6.6 


iVlCU 


iVlCU 


Mf»/1 




Q 

y 


7 
1 


2 


00Q07fi0901 


7.6 


+1.4 


39.3 


41.6 


41.4 


iVlCU 


iVlCU 


iVlCU 




SS 


22 




0145840201 


9.4 


-1.0 


19.7 


25.4 


25.7 


Thinl 


Thinl 


Thinl 




54 


38 


16 


0164561401 


10.0 


-0.2 


0.0 


32.0 


32.6 




Med 


Med 


0.92 


21 


13 


8 


0152832801+ 


10.3 


+0.3 


7.6 


9.4 


9.4 


Med 


Med 


Med 




13 


8 


5 


0152832901 t 


10.3 


-0.3 


5.6 


8.5 


8.5 


Med 


Med 


Med 




9 


6 


3 


0145840101 


10.4 


+0.1 


20.2 


26.2 


26.5 


Thinl 


Thinl 


Thinl 




33 


18 


15 


015283300lt 


10.5 


+0.0 


5.7 


8.2 


8.3 


Med 


Med 


Med 


0.93 


13 


10 


3 


0301270401 


10.7 


-0.2 


29.3 


39.0 


22.0 


Med 


Med 


Med 


0.96 


46 


32 


14 


0152833101+ 


10.7 


+0.3 


7.0 


8.7 


8.7 


Med 


Med 


Med 




9 


7 


2 


0024940201 


10.8 


+0.4 


22.5 


25.2 


25.2 


Med 


Med 


Med 


0.98 


50 


33 


17 


015283340lt 


11.0 


+0.3 


7.0 


8.7 


8.7 


Med 


Med 


Med 




22 


17 


5 


015283510lt 


11.5 


+0.0 


1.4 


4.8 


5.5 


Med 


Med 


Med 




3 


2 


1 


0152835501 t 


11.7 


-0.3 


1.7 


3.6 


3.5 


Med 


Med 


Med 




4 


3 


1 


01528357011" 


11.7 


+0.3 


8.0 


10.7 


10.7 


Med 


Med 


Med 




17 


10 


7 
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Table 1. Continued 



Observation 


Pointings 




Exposure ( 


ks) 




Filter 




FoV 




Sources 




ID 


I 


b 


pn 


MOSl 


MOS2 


pn 


MOSl 


MOS2 


Fraction 


total 


soft 


hard 


015283540lt 


11.8 


+0.0 


1.2 


3.0 


3.2 


Med 


Med 


Med 




3 


1 


2 


0150960301 


13.3 


— 1.3 


1.1 


4.8 


5.0 


Med 


Med 


Med 




3 


3 





0152835001+ 


14.8 


+0.0 


6.9 


8.5 


8.5 


Med 


Med 


Med 




20 


11 


9 


015283440lt 


15.0 


+0.3 


4.2 


8.7 


9.5 


Med 


Med 


Med 




11 


8 


3 


015283450lt 


15.0 


—0.3 


4.2 


8.9 


9.5 


Med 


Med 


Med 




11 


4 


7 


01 52834601 t 


15.2 


+0.0 


3.7 


6.1 


6.3 


Med 


Med 


Med 




7 


2 


5 


015283480lt 


15.3 


—0.3 


3.7 


7.0 


7.1 


Med 


Med 


Med 




8 


6 


2 


015283490lt 


15.5 


+0.0 


5.0 


8.4 


8.4 


Med 


Med 


Med 




18 


9 


9 


0040140201 


17.5 


—0.1 


12.3 


11.3 


12.3 


Med 


Med 


Med 




18 


7 


11 


0054540701 


18.0 


—0.7 


0.0 


23.1 


23.0 




Med 


Med 


0.98 


15 


6 


9 


0104460701 


19.2 


+0.0 


5.4 


8.4 


7.6 


Med 


Med 


Med 




11 


6 


5 


0051940101+ 


19.4 


+0.3 


4.8 


7.8 


7.8 


Med 


Med 


Med 




10 


5 


5 


0135745401+ 


19.6 


+0.0 


11.3 


13.7 


13.7 


Med 


Med 


Med 




13 


6 


7 


0051940301+ 


19.8 


+0.3 


6.1 


9.1 


9.1 


Med 


Med 


Med 




9 


3 


6 


0051940401+ 


19.8 


—0.3 


5.7 


9.1 


9.1 


Med 


Med 


Med 




16 


10 


6 


0104460301 


20.0 


+0.0 


7.7 


11.2 


11.0 


Med 


Med 


Med 




18 


10 


8 


0051940601+ 


20.2 


—0.3 


1.4 


3.8 


3.8 


Med 


Med 


Med 




2 


1 


1 


0051940501+ 


20.2 


+0.3 


6.8 


9.6 


9.6 


Med 


Med 


Med 




10 


4 


6 


0104460401 


20.4 


+0.0 


5.3 


10.3 


10.0 


Med 


Med 


Med 




9 


1 


8 


0135740701+ 


20.6 


+0.3 


4.8 


7.3 


7.4 


Med 


Med 


Med 




8 


2 


6 


0135745801+ 


20.6 


—0.3 


11.5 


13.2 


13.2 


Med 


Med 


Med 




15 


3 


12 


0104460901 


20.8 


+0.0 


11.3 


13.6 


13.6 


Med 


Med 


Med 




16 


3 


13 


0135740901+ 


21.0 


+0.3 


4.8 


7.8 


7.7 


Med 


Med 


Med 




11 


5 


6 


0135745901+ 


21.0 


—0.3 


5.8 


9.1 


9.5 


Med 


Med 


Med 




7 


2 


5 


0135745201+ 


21.2 


+0.0 


4.8 


7.1 


7.3 


Med 


Med 


Med 




7 


5 


2 


0122700801 


21.3 


—0.8 


12.5 


15.0 


14.8 


Med 


Med 


Med 


0.94 


14 


5 


9 


0112201001 


21.4 


+ 1.3 


0.0 


4.4 


4.5 




Med 


Med 




4 


2 


2 


0135745701+ 


21.4 


+0.3 


5.7 


7.4 


7.4 


Med 


Med 


Med 




10 


5 


5 


0135746801+ 


21.4 


—0.3 


2.7 


4.4 


4.4 


Med 


Med 


Med 




7 


2 


5 


0135746301+ 


21.6 


+0.0 


4.3 


6.7 


6.7 


Med 


Med 


Med 




6 


2 


4 


0135746401+ 


21.8 


+0.3 


1.7 


4.3 


4.5 


Med 


Med 


Med 




5 


1 


4 


0135741601+ 


21.8 


—0.3 


4.7 


7.7 


7.7 


Med 


Med 


Med 




11 


1 


10 


0135741701+ 


22.0 


+0.0 


3.6 


6.5 


6.8 


Med 


Med 


Med 




6 





6 


0135744801+ 


22.2 


+0.3 


1.1 


1.8 


2.0 


Med 


Med 


Med 




2 


1 


1 


0302560301 


23.2 


—0.3 


15.1 


17.9 


18.2 


Med 


Med 


Med 




22 


9 


13 


0400910101 


23.4 


+0.1 


10.2 


12.1 


12.1 


Med 


Med 


Med 




11 


3 


8 


0203850201 


24.3 


+0.1 


5.3 


6.8 


7.2 


Med 


Med 


Med 


0.89 


10 


6 


4 


0400910301 


25.4 


+0.1 


7.5 


9.1 


9.1 


Med 


Med 


Med 




20 


11 


9 


0301880401 


27.7 


+0.6 


8.4 


10.2 


10.1 


Med 


Med 


Med 




6 


3 


3 


0302970301 


28.4 


+0.1 


22.7 


26.7 


26.7 


Med 


Med 


Med 




28 


18 


10 


0301880901 


29.1 


+ 1.5 


14.9 


16.6 


16.6 


Med 


Med 


Med 


0.88 


17 


7 


10 


0046540201 


29.5 


+0.1 


6.7 


8.0 


8.4 


Med 


Med 


Med 




7 


3 


4 


0406140201 


29.8 


-0.7 


17.0 


18.5 


18.5 


Med 


Med 


Med 


0.98 


15 


4 


11 


0302970801 


30.3 


-0.5 


11.7 


15.5 


15.6 


Med 


Med 


Med 




28 


21 


7 


0203850101 


30.8 


+0.0 


20.6 


25.2 


25.2 


Med 


Med 


Med 


0.97 


31 


18 


13 


0203910201 


31.3 


-0.1 


5.4 


11.8 


12.2 


Med 


Med 


Med 


0.85 


6 


6 





0207010201 


31.1 


+0.2 


17.4 


18.7 


18.7 


Med 


Med 


Med 


0.89 


7 


1 


6 


0211080101 


32.4 


+2.1 


3.3 


6.8 


6.6 


Thinl 


Thinl 


Thinl 




7 


6 


1 


0136030101 


32.4 


+0.1 


17.8 


21.4 


21.9 


Med 


Med 


Med 


0.90 


27 


15 


12 


0306170201 


32.7 


+0.5 


9.9 


11.5 


11.5 


Thinl 


Med 


Med 


0.98 


7 


4 


3 


0017740601 


33.0 


+2.0 


20.0 


22.3 


22.3 


Med 


Med 


Med 




34 


15 


19 


0017740401 


33.1 


+0.0 


25.7 


27.7 


27.7 


Med 


Med 


Med 


0.77 


16 


8 


8 


0302970201 


36.7 


-0.1 


2.7 


10.0 


13.3 


Med 


Med 


Med 




6 


3 


3 


0136030201 


38.5 


+0.0 


13.3 


16.0 


16.0 


Med 


Med 


Med 


0.98 


10 


6 


4 


0305580201 


43.0 


+0.8 


22.3 


27.8 


28.2 


Med 


Med 


Med 


0.89 


22 


6 


16 



+ Observations corresponding to the Galactic X-ray Survey Programme (XGPS) 



6 R.S. Warwick, D. Perez-Ramirez and K. Byckling 



hard bands were obtained as the straight sum of the rates recorded 
in the respective input bands. We also calculated an effective detec- 
tion maximum likelihood (maxl) in the soft and hard bands by the 
prescription given in the SAS emldetect documentation (see also 
iMateos et alTlOOSl) . 

A further compression of the information for each detected 
source was achieved by combining the count rates measured in the 
different cameras for the soft band and, separately, for the hard 
band. For the count-rate measurements for a particular source in 
a given camera to be considered valid, a requirement was that the 
2XMMi source parameter, frac, should be greater than 0.8 for that 
camera {i.e., that at least 80% of the source response was contained 
within the camera's active field-of-view). Count-rate measurements 
passing this criterion were included irrespective of whether the 
source was classed as detected or not in the camera/band combi- 
nation under consideration (see below). The actual compression 
process involved first averaging valid MOS-1 and MOS-2 count 
rates using a statistical weighting to give a combined-MOS count 
rate estimate and then scaling the latter measurement to pn equiv- 
alent unitfl Finally the scaled-MOS and pn measurements were 
averaged, again using a statistical weighting, to give the final soft- 
and hard-band count rates. In following this process, the resultant 
source count rates may derive from one, two or three measurements 
depending on which combination of cameras were producing valid 
data for the source. Hereafter in this paper we quote count rates in 
the soft, hard and total (soft+hard) bands for the combined EPIC 
cameras in units of pn count ks"^. 

The first step in excising relatively low-significance sources 
from the full 2XMMi source list involved the setting of a threshold 
value of maxl = 10 for the soft/hard band detection. Sources were 
only retained in the source list if this threshold was exceeded in 
either the soft or hard bands (or both) and in one or more of the 
EPIC cameras. For a source detection to be classed as valid in a 
particular camera, we required frac > 0.8 in that camera. Further 
criteria for source removal included the 2XMMi parameter Extent 
> {i.e., the exclusion of all non-pointlike sources) and the SUM- 
FLAG > 3 (i.e., a quality check to exclude sources described as 
"located in an area where spurious detections may occur and may 
possibly be spurious" - see lWatson et al.ll2009l) . 

After filtering the original 2XMMi source list using the above 
methodology, we visually inspected the soft- and hard-band images 
of each field (produced from the standard data products available 
from the XMM-Newton archive). As a result a number of sources 
classified at this stage as "valid detections" were flagged for re- 
moval. Frequently such source exclusions stemmed from the con- 
fusion of one source with another, generally brighter, object. A few 
sources embedded in bright diffuse regions were also removed. At 
this point we also excluded regions encompassing an obvious bright 
target source - essentially by spatially masking areas within the ob- 
servation field of view. On occasion, regions affected by bright ex- 
tended objects or by the scattering rings from very bright out-of- 
field-of-view sources were also excluded. Observations in which 
the target region or other specific regions have been excluded are 
flagged in Table I. 

Implementation of the above selection process left us with a 
catalogue of ~ 2700 sources. A source could be included in this 
list based on a nominal detection (above a maximum likelihood 



The scale factor for converting the combined-MOS to pn count rates was 
determined for each band on the basis of the average pn:combined-MOS 
count rate ratio derived for the full source sample 



Table 2. Source detection statistics. 



Camera 


Soft-Band 
Only 


Hard-Band 
Only 


Both Soft 
and Hard 


All 


pn only" 


349 


222 


18 


589 


MOS-1 only" 


55 


57 


3 


115 


MOS-2 only" 


56 


67 


7 


130 


Multiple'' 


626 


453 


291 


1370 


Total 


1086 


799 


319 


2204 



" Detections in a single camera only 
Detections in at least two cameras 



threshold of 10) in just one band and one camera. Visual inspection 
suggested that at the limit, such sources were typically quite faint 
and in some cases potentially unreliable. Since a detailed analysis 
of the reliability of 2XMMi sources as a function of the detection 
maximum likel ihood threshold wa s beyond the scope o f the investi- 
gation (but see'Watson et a l.ll2009l : lMateos et al.ll2009h . we instead 
employed an empirical approach to arrive at a robust final source 
catalogue based on a further signal-to-noise selection. More specif- 
ically, we calculated the ratio of the count rate to count-rate error 
for each source using the combined camera measurements. We then 
set a minimum requirement that this ratio should exceed a value of 
4 in the soft and/or hard bands. This further selection resulted in the 
exclusion of 252 sources detected only in the pn camera and 101 
(133) sources seen only in the MOS-1 (MOS-2) cameras, with just 
10 sources removed for which there was a simultaneous detection 
in more than one camera. 

After applying all the above, we were left with a catalogue 
of 2204 point sources drawn from the 116 survey fields. Table |2] 
summarises the source detection statistics in terms of the soft/hard 
bands and the pn/MOS-l/MOS-2 cameras. 

Although the field selection described in §2 involved the ex- 
clusion of duplicate observations of the same target or pairs of ob- 
servations with fields-of-view which significantly intersect, some 
modest degree of field overlap was permitted. As a consequence of 
this overlap, there are a total 44 source duplications in the total of 
2204, i.e., 2% of the sample. Within these 44 sources, 30 of which 
are soft detections and 14 hard, there are two instances of the same 
source being twice duplicated. In the following analysis we have 
ignored this duplication and treated all 2204 entries as individual 
sources. 



4 PROPERTIES OF THE SOURCE SAMPLE 
4.1 Spatial distribution 

The distribution in Galactic longitude and latitude of the sources 
which comprise our serendipitous source sample is illustrated in 
Fig [3] These are the observed distributions without any correction 
applied for sky coverage or for the depth of the underlying obser- 
vations. 

It is evident from Figs. [T] and [3] that the coverage of the 
Galactic plane by XMM-Newton is far from uniform. Near the 
Galactic Centre, within |i|<2°, there is an obvious concentration 
of pointings. Similarly the regions near 1 = 7°, 11°, 15° and be- 
tween 19°-22° have relative good coverage in narrow strips along 
the Plane. The l atter c orresponds to the "XGPS region" studied 
by iHands et all (l20oi). whereas the other three directions have 
enhanced coverage as a result of later phases of the XGPS pro- 
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Galactic Longitude 



Galactic Latitude 



Figure 3. a) The source distribution in Galactic longitude. No corrections 
have been applied for the number or depth of the observations in the survey; 
b) The same in Galactic latitude. 



gramme. The observations drawn from the XGPS programme are 
flagged as such in Table [T]- generally these consist of sequences 
of relatively short exposures with some overlapping of the fields of 
view so as to give extended, albeit shallow, sky coverage. 



4.2 Flux calibration and distribution 

The distribution of measured total count rate for the source sam- 
ple is illustrated in Fig. |4] Although the full distribution extends 
over 3 decades, 94% of the sources have total count rates in the 
range 2-100 pn ct ks~^. Figure [5] shows how the total count rate 
splits between the soft and hard bands. Roughly 50% of the sources 
are detected only in the soft band, with 36% detected only in the 
hard band and just 14% detected in both bands (see Table O. 
With the objective of dividing the sample into non-overlapping 
soft and hard subsets, we calculate a broad-band hardness ratio, 
HR = (H - S)/(H + S), where S and H refer to the soft-band 
and hard-band count rates respectively. We then use HR=0, corre- 
sponding to the diagonal line in Fig. |5] (i.e., equal soft and hard 
count rates), as the boundary between the two designations. This 
provides a close match to the soft-only and hard-only detections; 
on the other hand, for the sources detected in both bands, the di- 
vision somewhat favours the hard category. Hereafter we refer to 
sources as soft or hard depending on their position relative to the 
HR = fiducial. On this basis there are 1227 soft sources and 977 
hard sources in our sample. 

Given the very different characteristics of the underlying 
source populations (see below), different spectral forms have been 
assumed for the soft and hard sources in deriving factors to convert 
the measured count rates to energy flux. In both bands we apply a 
spectral model consisting of a power-law continuum of photon in- 
dex r subject to line-of-sight absorption in a column density. Nr. 
For the soft sources, the spectral parameters were set to F = 2.5 and 
= 10^^ cm~'^ resulting in an energy conversion factor (ECF) of 
4.9x10" ct cm'^ erg ^ , where both the count rate and unabsorbed 
energy flux relate to the 0.5-2 keV band. For the hard sources, we 
use r = 1.0 and A^H = 3 X 10^^ cm'^ leading to an ECF of 7.9x lO"' 
ct cm^ erg^^, applicable to the 2-10 keV band. For a source with 
a nominal count rate of 2 pn ct ks~^ in the soft band, the equiva- 
lent unabsorbed flux is 4x lO"^'^ erg s"^ cm"^ (0.5-2 keV). The 
same count rate in the hard band equates to an unabsorbed flux of 



o 
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Figure 4. The distribution of the total (soft+hard) count rate for the source 
sample in units of pn ct ks^^ (0.5-12 keV). 



in O 
J£ O 



c 
a. 



TTTTl 1 1 I I lllll 1 1 I I I nil 1 r-TTTTTTT— 1 1 I I I 




2.5 x 10"'* erg s"" cm"^ (2-10 keV). 



Soft Count Rate (pn ct/ks) 

Figure 5. Soft-band count rate versus hard-band count rate for the source 
sample. Soft-band only detections are shown as red crosses and hard-band 
only detections as blue crosses. Sources detected in both the soft and hard 
bands are shown as black circles. The thick dashed line con'esponds to a 
hardness ratio HR = (see text). Count rates below 0.1 pn ct/ks are plotted 
at this value. 



4.3 X-ray Spectral properties 

Fig.[6]shows the distribution of the broad-band hardness ratio, HR, 
for the source sample. The distribution has a saddle-like form, with 
a broad minimum centred around HR=0 separating maxima at, or 
near, the two extremes. We note that our focus on source detections 
in either the soft or hard bands (or both) might introduce a bias 
against sources of intermediate hardness in which the count rate 
is spread fairly evenly across the two spectral channels. However, 
when we plot the distribution in Fig. |6] including only relatively 
bright sources (e.g., sources with a total count rate greater than 
20 ct/ks) the saddle distribution remains broadly unchanged. We 
conclude that the underlying source population splits rather cleanly 
into the designated soft and hard categories. 

We have also investigated the spectral distribution of the 
sources with respect to the four 2XMMi bands, namely Bands 2-5. 
In this context a standard approach is to use X-ray two-colour dia- 
grams, where each axis represents a different hardness ratio calcu- 
lated for a pair of (adjacent) bands. However, for our sample which 
is strongly affected by absorption, the tendency was for such dia- 
grams to be either dominated by errors (due to low intrinsic count 
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-0.5 



Broad— band Hardness Ratio 

Figure 6. Distribution of tlie sources witli respect to tlie broad-band hard- 
ness ratio, HR. 



rates in both selected bands) or, at least, very compressed along 
one axis. We have avoided this problem by calculating instead a 
set of Band Index (BI) values, where BJ-N is calculated simply 
as the count rate in Band N normalised to the total source count 
rate (i.e., in this case the sum of the count rates in Bands 2-5). 
The comparison of BI values for adjacent bands - see Fig. |7] then 
provides very similar information to a set of two-colour hardness 
ratio diagrams. W e note that "quantile" analysis, as described by 
iHong et ai] ( |2004|) . represents an alternative strategy for investigat- 
ing the spectral properties of source samples in the low count rate 
limit. 

Each BI diagram in Fig.|7]shows a smoothly varying distribu- 
tion representing the range of source characteristics from the hard- 
est through to the softest spectral forms. A rough calibration of this 
spectral diversity is provided by the model curves in Fig.|7] which 
correspond to power-law spectra with F ranging from 0.5-3, sub- 



ject to absorption ranging from Nh 



10^ 



cm . However, 



the sources with the softest spectra are not quite accommodated by 
the limits T = 3 and Nh = 10^° cm"^ in the BI-2 versus BI-3 
diagram and are, in fact, better represented by a thermal (Mekal) 
spectrum with the same limiting Nh and kT ~ 0.5 keV (the 
dashed curve in the top panel of Fig.|7j. From the BI-3 versus BI- 
4 diagram, it is evident that transition from our soft to hard charac- 



terisation occurs for Nh = 3 x 10 



(depending 



on the slope of the continuum). In the hard spectral regime repre- 
sented by the BI-4 versus BI-5 diagram, one finds that the bulk of 
the signal is contained in Bands 4 and 5 for Nh > 3 x 10^^ cm~^. 
Evidently, there are a significant number of sources with BI-5 
> 0.7, implying Nh in excess of 10^^ cm~^, provided the un- 
derlying continuum is not exceedingly flat (F ^ 0.5). 

4.4 Column density constraints 

We have used the measured BI values to set more specific con- 
straints on the line-of-sight column density Nh for each source. 
Assuming an absorbed power-law spectral model, we calculated a 
two-dimensional array of BI values (one such array for each of 
the four spectral bands 2-5), with log A^'^f along the x-axis (log 
Nh = 20 — 24 cm~^ in steps of 0.1) and the assumed F along 
the y-axis (F — 0.5 — 2.5 in steps of 0.1). For a given source, we 
identified the region in each array bounded by the measured BI 
value and its error bar (we actually used 1.65a errors correspond- 
ing roughly to 90% confidence combined in quadrature with a 2% 
systematic). The full constraint was then derived by taking the min- 
imum bound on log Nh consistent with all four BI plots. We set 
the final estimate for log Nh as the mid-point of the bound and 




0.5 

BAND INDEX 2 




0.5 

BAND INDEX 3 




0.5 

BAND INDEX 4 

Figure 7. X-ray two colour diagrams in the form of Band Index (BI) plots. 
Here Band Index N is simply the count rate in band N divided by the total 
count rate summed across the 2XMMi Bands 2-5. Sources categorised as 
soft and hard are shown as red and blue respectively. For clarity only sources 
with a signal-to-noise threshold in excess of six are plotted. The solid curves 
represent the locus for a given photon index, F, as the column density Nh 
varies from 10'^'' cm~^ to 10^* cm~^ (with the dots marking in an anti- 
clockwise order: 10^^ 3 X lO^^; 10^^; 3 X 10^^ and 10^^ cm-^). The 
values of F are as indicated (ordered such that the higher F loci are at a 
further distance from the origin). The top panel also shows the locus (dashed 
curve) for a thermal Mekal component with kT = 0.5 keV, with the same 
Nh range as above. The dashed diagonal lines represent an upper bound at 
which the full count rate is accounted for by the two bands in question. 
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took the half-width of the bound as the error. This simple recipe 
worked weljf] when applied to the hard source sample, but was of 
more limited value for soft sources, where the range of power-law 
models represents a poor match to the underlying (thermal) spectra. 

The result of applying the above procedure to the hard sam- 
ple is illustrated in Fig. [8ja), which shows the X-ray derived Nh 
values versus the source broad-band hardness ratio, HR. Although 
the errors for individual sources are typically quite large (±0.45 
in log Nh, on average), there is a clear trend with Nh ranging 
from ^ 10^^ cm^^ to > 10^"^ cm^"^ as HR varies from 0-1. Fig. 
[Stb) shows the comparison of the column density derived from the 
X-ray measurements with the total Galactic column density in the 
source direction, as inferred from the inf rared dust m easurements 
of COBE/DIRBE and IRAS/ISSA (Schlegel et alJl99S) . From this 
plot it is evident that, for hardness ratio values HR ^ 0.8, which 
represents 50% of the hard sample, the X-ray Nh is very compa- 
rable to the total Galactic column density. In this regime, the X-ray 
source could be either embedded in a localised high density region 
(perhaps intrinsic to the X-ray source), located at a significant dis- 
tance across the Galactic Plane or be an extragalactic interloper. We 
consider the latter possibility in section l6^ 

Finally, we note that the X-ray Nh values derived for the hard 
source sample show a clear dependence on Galactic longitude and 
latitude, albeit with a large point-to-point scatter. Sources located 
within 2.5°of the Galactic Centre have Nh values typically a fac- 
tor 2 higher than those found outside of this region, the average 
in the Galactic Centre direction being Nh = lO^'^ * cm~^. Simi- 
larly, sources located within l°of the Galactic plane typically have 
four times the column density of sources which are 1°- 2.5°off the 
plane. 



5 NEAR-INFRARED COUNTERPARTS 

The identification and characterisation of longer- wavelength coun- 
terparts represents a crucial step in determining the nature of Galac- 
tic X-ray sources. The ideal starting point for this would be a sub- 
arcsecond imaging survey in a set of wavebands not too strongly in- 
fluenced by Galactic absorption, with commensurate high precision 
astrometry. At present a comprehensive dedicated image database 
of this quality is not available for the XMM-Newton survey fields. 
However, we have conducted a systematic investigation based on 
cross-correlations with two near-infrared surveys which are avail- 
able, namely 2MASS and UKIDSS. 



5.1 Cross-correlation with 2MASS 

The Two Micron A ll Sky Survey (2MASS) dCutri et all I2OO3I : 
ISkrutskie et al.ll2006[) provides uniform coverage of the entire sky 
in three near-infrared (NIR) bands, namely in J (1.25/i), H (1.65/i) 
and Ks (2.17/^) on a scale of 2.0"(this is the pixel size of 2MASS 
images; the FWHM of the PSF was typically 2.5"-3.0")- At high 
latitude the nominal survey completeness is 15.8, 15.1 and 14.3, 
respectively. However, in the Galactic Plane the measured source 
counts turn down at limits one magnitude (or more) brighter, be- 
cause of the effects of source confusion on the detection thresholds. 



' The procedure is, of course, only an approximation to detailed model 
fitting of the available X-ray spectral data. Nevertheless it provides very 
useful indicative results. 




Broad— band Hardness Ratio 

Figure 8. a) The line-of-sight Nh estimated from the BI values plotted 
against the broad-band hardness ratio, HR. Only sources categorised as hard 
and with a signal-to-noise threshold > 6 are shown. An underlying power- 
law continuum spectrum, with F in the range 0.5-2.5, is assumed, b) A 
comparison of the X-ray derived Nh with the total Galactic Nh in the 
source direction determined from infrared dust measurements. The y-axis 
represents the difference between the X-ray and the infrared estimates of 
log Nh- 
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X-ray Position Error (arcsec) 

Figure 9. The distribution of X-ray (Icr) position errors for the source sam- 
ple. 

The astrometric uncertainty is generally less than 0.2", although 
this may be compromised in very confused regions. 

As a first step, we cross-correlated our X-ray source sample 
with the 2MASS catalogue and, for each X-ray source position, ex- 
tracted the set of 2MASS sources within a radius of 20". A total of 
27485 2MASS sources were associated with the 2204 X-ray source 
positions via this process. Here we use the X-ray source position 
and position error quoted in the 2XMMi catalogue for the indi- 
vidual source observations (parameters: RA, DEC, poserr), where 
the position error includes a systematic error added in quadrature 
with the statistical error (the former was set to 0.35" if the field 
astrometry was corrected by reference to the USNO Bl.O optical 
catalogue, or l"otherwise). Fig|9] shows the distribution of X-ray 
position errors for the source sample; the range is from 0.37"to 3", 
with an area- weighted average of 1.5". In the event we selected 
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X— ray Error Circle Radius 

Figure 10. The cumulative probability of finding an associated 2MASS 
source in an X-ray error circle, plotted as a function of the error circle 
radius (in units of the la X-ray position error). The analysis is restricted 
to 2MASS sources with J = 5 — 15 found in the fields of the soft X-ray 
sources. For the current work, the radius of the X-ray error circle was taken 
to be 2.5 times the Icr X-ray position error (as indicated by the vertical dot- 
ted line), thereby encompassing around 90% of the likely NIR counterparts. 

only those X-ray sources with X-ray position errors <2"for the 
cross-correlation study, resulting in a reduced sample size of 2016 
X-ray sources. 

A preliminary investigation of the incidence of 2MASS 
sources at or near the X-ray positions demonstrated that, at least for 
the soft X-ray source sample, there was a significant excess over the 
number expected by chance, presumably reflecting the presence of 
a real counterpart in a substantial number of casefl The distribu- 
tion of this excess population with respect to the X-ray source po- 
sitions was studied by considering how the net number of sources 
contained within the soft-source error circles (over and above the 
background rate) varied as a function of the assumed X-ray error 
circle radius (FigllOb. On the basis of this empirical information, 
we set the X-ray error circle radius to be 2.5 times the X-ray posi- 
tion error quoted in the 2XMMi catalogue, at which point 90% of 
potential real counterparts will be contained within the error circle. 

For the soft X-ray sources, the number of error circles contain- 
ing a 2MASS source was 987 out of 1120 (88%), whereas for the 
spectrally hard sources the statistics were 393 out of 896 (44%). 
In most instances, there was only one 2MASS source within the 
X-ray error circle, but for the 85 X-ray sources for which this was 
not the case, we treated the brightest 2MASS source (in J band) 
as the potential counterpart. Figure [TT] shows how the incidence 
rate of 2MASS stars within the X-ray error circles accumulates as 
a function of the 2MASS J magnitude for both the soft and hard 
source categories. The chance coincidence rate, calculated using 
the 2MASS catalogue information for the field regions within 8~ 
20" of the X-ray source position, is also indicated. All three curves 
begin to flatten and eventually reach a plateau, around J ~ 16, con- 
sistent with the fact that the spatial density of sources per unit mag- 
nitude in this set of Galactic fields peaks at J ~ 15 (i.e., at a level 
roughly a magnitude brighter than the completeness limit of the 

* In establishing the coincidence rate of 2MASS sources with the soft and 
hard X-ray source samples, we did not employ any discrimination based 
on the 2MASS quality flags. However, the flags were use to select sources 
with good photometry in the later investigation of the impact of interstellar 
reddening on the stellar colours. 
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Figure 11. The cumulative fraction of X-ray error circles containing a 
2MASS star as a function of the J magnitude of the star. Separate curves 
are shown for the soft and hard X-ray source samples. The chance coinci- 
dence rate for field stars is shown by the filled histogram. 
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Figure 12. The fraction of error circles containing 2MASS stars as a func- 
tion of the J magnitude for the soft (top panel) and hard (bottom panel) 
source samples. This is a differential plot, i.e., the ordinate is the number 
per magnitude bin normalised to the total number of sources in the soft/hard 
sample. The predicted distribution of chance coincidences with field stars, 
based on the statistics for 2MASS stars located at radial offsets between 
8"-20" from the X-ray position, is shown as the filled histogram. This 
prediction takes into account the declining fraction of error circles which 
are available at faint magnitudes due to the presence of brighter stars. 



2MASS survey at high latitude, as noted earlier). The plateau level 
for soft and hard distributions is as quoted above, whereas for the 
chance distribution it is ~ 30%. 

Figure [12] shows the derived magnitude distribution for the 
brightest 2MASS stars found to be positionally coincident with 
the X-ray source positions. The predicted contamination due to 
chance coincidences is also shown, where a correction has been 
incorporated for the fact that an error circle may not be "avail- 
able" at faint magnitudes if a brighter star is already present. For 
the soft X-ray source sample, the magnitude distribution extends 
from J=7-16. Down to J=16, the underlying contamination due to 
field stars amounts to roughly 9% out of the 88% of error-circles 
which contain 2MASS stars down to this magnitude limit. In con- 
trast, the magnitude distribution for the hard sample shows a much 
more modest excess of coincidences with respect to the predicted 
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Broad— band Hardness Ratio 

Figure 13. The fraction of sources with a 2MASS source with J< 16 in 
the eiTor-circle as a function of HR. The chance coincidence rate is also 
indicated. 



chance rate; for J < 16, the coincidence rate is 38% of which 24%, 
i.e., roughly two-thirds, can be attributed to contamination by field 
stars. 

The probability of finding a positionally coincident 2MASS 
source is clearly very different for the soft and hard X-ray source 
samples. We have investigated this in more detail by considering 
how this probability varies as a function of the broad-band hard- 
ness ratio, HR. Fig. [T3] shows the result. The probability of finding 
a 2MASS star brighter than J = 16 in the X-ray error circle (inclu- 
sive of chance hits) exhibits a step function behaviour, remaining 
at close to 90% for the spectrally soft sources with HR in the range 
- 1 to -0. 1 , followed by a rapid transition to a level of about 50% for 
HR values near followed by further decline to a level very near to 
the chance rate for the sources with the hardest spectra. 



5.2 Cross-Correlation with UKIDSS 

For our application, the 2MASS survey becomes incomplete below 
J ~ 15. In contrast the United Kingdom Infrared Deep Sky Sur- 
vey (UKIDSS, Lawrence et al. 2007) provides better spatial reso- 
lution (0.2" pixel scale with typically ^ l"seeing) and is signifi- 
cantly deeper. The nominal limiting magnitudes estimated for the 
UKIDSS survey are 1=19.77, H=19.00 and K=18.05, with uncer- 
tainties of up to ~ 0.2 magnitudes. As with 2MASS, the appro- 
priate limits for the Galactic plane will be considera bly brighter 
depen ding on the actual crowding within a given field jLucas et ahl 
l2008b . The photometric data used in this analysis is taken from the 
UKIDSS Galactic Plane Survey (GPS) Data Release 4. 

UKIDSS GPS I-band data were available for a total of 1 1 17 X- 
ray sources drawn from our sample of X-ray sources with position 
errors less than 2", with the coverage split roughly equally between 
the soft and hard samples {i.e., 596 soft and 577 hard). For this 
analysis we extracted a preliminary set of UKIDSS stars within a 
nominal 10" of each X-ray position, leading to a preliminary set of 
36974 stars down to J=20. Using the same X-ray error circle radii as 
for the 2MASS study, we repeated the analysis of the coincidence 
rates of UKIDSS stars within the X-ray error circles. As before, we 
make the assumption that the brightest star in the circle (in J band) 
is the most likely counterpart. 

As we push down to a J magnitude of 20, the incidence rate 
of infrared stars in the X-ray error circles rises to 97% for the soft 
sources and 92% for the hard sources. Of course at these faint lev- 
els, given the stellar density and size of the error circles, the chance 
rate is commensurately high, namely 85%. Figure[T4]shows the re- 
sulting magnitude distributions for the brightest UKIDSS star in 
the X-ray error circle. 

For the soft sources, the magnitude distribution is truncated 
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Figure 14. The fraction of error circles containing UKIDSS stars as a func- 
tion of the J magnitude for the soft {top panel) and hard {bottom panel) 
source samples. This is a differential plot, i.e., the ordinate is the number 
per magnitude bin, normalised to the total number of sources in the soft/hard 
sample. The predicted distribution of chance coincidences with lield stars, 
based on the statistics for UKIDSS stars located at radial offsets between 
7.5"-10" from the X-ray position, is shown as the filled histogram. This 
prediction takes into account the declining fraction of error circles which 
are available at faint magnitudes due to the presence of brighter stars. 



below J < 10 as a result of saturation effects in the UKIDSS 
images, with the apparent sub-peak near J « 11 presumably at- 
tributable to the same effect. As was evident in the equivalent 
2MASS distribution, there is a sharp cutoff between 1=15-16. De- 
spite the fact that the UKIDSS survey goes between 3-4 magni- 
tudes deeper than 2MASS, the evidence is that very few real soft- 
source counterparts have J > 15.5. For the soft sources, the coin- 
cidence rate with bright UKIDSS stars (J < 16) is 91% inclusive 
of a chance rate of 9%, values which are fully consistent with the 
equivalent 2MASS estimates. 

For the hard sources, the coincidence rate with bright UKIDSS 
stars (J < 16) is 45%, which splits into a likely chance rate of 29% 
and an excess rate {i.e., potential real counterparts) of 16%. These 
percentages are again very comparable to those obtained from the 
2MASS cross-correlation. There is also evidence that real counter- 
parts are seen at fainter magnitudes, with the J=I6-20 coincidence 
rate of 47% comprising a likely chance rate of 38% and an excess 
rate of 9%. It follows that ~ 25% of the hard source population 
may have counterparts visible, as the brightest star in the error cir- 
cle, down to J=20. Of course, for every error circle for which this 
is the case, there are roughly 2.7 instances where the brightest star 
is just a chance positional coincidence of a field star. Clearly this 
statistic implies detailed follow-up of the hard-source population 
carries a large overhead. 



5.3 Reddening of potential stellar counterparts 

In order to investigated the colours of the NIR stars found inside the 
X-ray error circles, we have carried out a further selection aimed at 
removing those objects affected by known issues in relation to their 
photometry. The 2MASS catalogue includes four quality flags. For 
our purpose, we utilise the so-called the photometric QUALITY 
flag and select only those stars flagged as AAA {i.e., valid mea- 
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Figure 15. Top panels: The H-K versus J-H colours of the brightest stars 
in the X-ray errors circles, shown separately for the soft and hard X-ray 
source samples. The curved line is the locus of dwarf stars from FOV to 
M6V. The diagonal line shows the reddening vector for Ay = 20. Bottom 
panels: The J-K colour versus J magnitude of the brightest stars in the X-ray 
errors circles, shown separately for the soft and hard X-ray source samples. 
The diagonal line shows the impact of Ay = 20 on the J-K colour and J 
magnitude of a given star. 




AV 

Figure 16. A comparison of the visual absorption Ay determined from 
the 2MASS/UKIDSS stellar colours with the estimated line-of-sight gas 
column density Nh derived from the X-ray spectral characteristics. The 
curve shows the standard relationship Nh = 1-79 X 10^^ Ay cm^^ 
jPredehl & SchmitJI 19951) . 



surements with signal-to-noise ratio >10 in J, H and K). Similarly 
for the UKIDSS stars we required the error bit flag, ppErrbits < 
256, for each of the three near-IR bands and also pstar ^ 0.99 
(i.e., the removal of sources with non-stellar profiles). A total of 
840 2MASS stars, representing the brightest object in the X-ray er- 
ror circle, passed this selection step, with the equivalent number for 
UKIDSS being 444. A total of 130 stars were common to both lists 
and, in the event, we used only the 2MASS data for this subset of 
sources. 

Figure [73] shows the H-K versus J-H two-colour diagram for 
the combined 2MASS/UK[DS^| sample of coincident stars, with 



the results for the soft and hard X-ray sources shown se parately. 
The l ocus of main- sequence dwarf stars from FOV to M6V jPicklej 
Il998l) and the reddening vector for a visual absorption Ay = 20 
are also defined. For the stars potentially associated with soft X- 
ray sources, there is a dense grouping in the region of the two- 
colour diagram encompassed by the locus of non-reddened main 
sequence dwarf stars, with a spread along the reddening direction 
largely encompassed by an Ay of up to 20. In the case of the hard 
X-ray sources, the implied reddening is much more substantial for 
many of the stellar candidates, although the bulk of the popula- 
tion are still bounded by Ay < 20. Particularly for the hard band 
sample, there are a number of stars which lie off to the right of 
the main distribution (i.e., have very red H-K colours); these may 
be pre-main sequence (T Tauri) stars exhibiting an infrared excess 
(rf lKoenigetal]|2008l) . 

Figure [TS] also shows the variation of the J magnitude versus 
the J-K colour, with the impact of a visual absorption Ay ~ 20 
again indicated. The absolute magnitude, Mj, for m ain sequence 
stars of type FOV to MOV ranges from 2.43-5.72 jCovev et al] 
l2007l) . i.e., a A J ~ 3.3 for unreddened stars at a fixed distance. If 
we compare this with the spread of A J « 7 apparent in the colour- 
magnitude diagram for the bulk of the stars associated with the soft 
X-ray sources, then the implied scatter in distance is roughly a fac- 
tor of 5. Of course, the presence of significant numbers of dMe stars 
of spectral class later than MO would compromise this argument. 
The stars found in the hard X-ray source error circles are signif- 
icantly fainter and redder than the population associated with the 
soft X-ray sources. However, as noted earlier, a substantial fraction 
of the stars linked to the hard X-ray sources will be chance coinci- 
dences with field stars and to first order the scatter in the both the 
H-K versus J-H two-colour diagram and in the J-K versus J colour- 
magnitude diagram for the hard sources must reflect the underlying 
properties of the field star distribution. 

For those stars linked to hard X-ray sources, we have esti- 
mated the visual absorption by projecting the star's position in 
the two colour diagram onto the reddening vector and determining 
the offset relative to an origin at J-H=0.61, H-K=0.11 (the colours 
of an unreddened GV star). In Fig. [T6l the resulting estimate of 
Ay is plotted versus the X-ray column density derived in §4.4. 
Compared to the standar d conversion, Nh = 1.79 X 10^^ Ay 
jPredehl&Schmittll 19951) . most of the points in the diagram lie at 
Nh values significantly greater than that implied by the Ay de- 
termination. This is to be expected if many of the associated NIR 
objects are, in fact, foreground stars, albeit sufficiently distant to 
be subject to significant reddening. One way of picking out poten- 
tially true counterparts might be to require the derived Ay and A^j^ 
estimates to be comparable; unfortunately the large uncertainties 
implicit in the estimation of both parameters mitigates against this 
as a practical scheme. 



6 DISCUSSION 

6.1 Nature of the Soft X-ray Source Population 

In the soft band, the nominal sensitivity limit of our survey is 
around 2 pn ct/ks corresponding to an unabsorbed flux. Fx = 4 x 



10" 



erg cm s (0.5-2 keV). This is roughly 5 times deeper 



^ We have not ap plied the photometric transformations detailed in 
iHodgkin et alU200(il) to convert the UKIDSS measurements to the 2MASS 



photometric scale, since the scatter in the plots far exceeds the magnitude 
of the corrections. 
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than the ROSAT survey reported bv lMorlev et"al] j200lh. which in 
turn was a factor 5 deeper than the Einstein GPS lliertz & Grindlavl 
1984; Hertz & Grindlav '1988') and the ROSAT GPS (^Motch et al ' 
1991; Motch et al. 1997). By way of comparison, the stellar sur- 
vey carried out recently under the auspices of th e Extended Chan- 
dra Multiwavelength Project jCovev et al.ll2008i) reaches similar, 
or somewhat deeper, limits to those reported here. These surveys 
and many other studies have established that the Galactic source 
population in the soft X-ray band largely comprises of late-type 
stars with active coronae. More specifically the ROSAT all-sky sur- 
vey clearly demonstrated that main-sequence stars from early-F to 
late-M dominate the source statistics and that stars younger than 
~ 2 Gyr are typically more luminous in X-rays than older popula- 
tions; this is a consequence of the reducing efficiency with which 
the corona is heated by the dy namo mechan ism as the rate of stellar 
rotation declines with age (see lGudelll2004l and references therein). 
At high galactic latitude, as the flux limit is lowered, the surveyed 
volume will eventually extend beyond the scale height of the young 
star population, thereby altering the balance betwee n relatively old 
and relatively young stars (e.g., iMicela et al.]|2007h . Although this 
is not a consideration which applies directly to a narrow Galactic 
plane survey, line-of-sight absorption will suppress the surveyed 
volume most severely for higher Lx sources, presumably giving 
rise to a similar bias. 

From the cross-correlation of our soft X-ray source sample 
with the 2MASS catalogue, we concluded that 88% of the error 
circles contain stars with J < 16, with only 9% of this figure at- 
tributable to chance coincidences. Furthermore, since our error cir- 
cles were scaled so as to contain ~90% of the true counterparts, 
it follows that active stellar coronae must account for almost our 
entire sample of soft X-ray sources. When the search for counter- 
parts was extended some 4 magnitudes fainter through the cross- 
correlation with the UKIDSS catalogue, the evidence for a cut-off 
at J f»16 was quite starljfl This prompts the question of why the 
stellar identifications are restricted to bright magnitudes? 

In Fig. [it] we explore the distances out to which main se- 
quence stars of different spectral class may be detected in both 
the NIR and our XMM-Newton survey. In the NIR, the constraint 
J ^ 16 sets a distance limit which varies from 2.7 kpc for early F 
stars down to ~ 70 pc for late M stars. In contrast if we assume 
stellar X-ray luminosities typical of normal solar-type sta rs, as rep- 
resented by the NEXXUS survey jSchmitt & Liefkei2004h . then the 
X-ray survey horizon is an order of magnitude closer. Of course, as 
the level of coronal activity increases, so does the surveyed volume 
in the X-ray band. As is clear from Fig.[T7] if the X-ray horizon is 
to match the J=16 boundary, the X-ray luminosity will need to have 
increased to the point at which stellar-coronal emission is known to 
saturate, namely at Lx ~ x Li,oi (Gudel 2004 and references 
therein). This, almost certainly, explains the origin of J ~ 16 cut-off 
alluded to earlier. 

It is true that previous X-ray surveys have identified some 
coronal sources emitting above the nominal saturation limit, 
one explanati on being that these are systems cau g ht whilst flar- 
ing (e.f>., Fleming e t aL||l995l ; iMorlev et alJboOll) . 'Fleming et al. 
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1995,) estimate that ~25% of solar-like stars and -^50% of dMe 

^ If we compare the corresponding X-ray and infrared flux limits. 
Fx = 4 X 10~^^ erg cm~^ and J = 16, through the relation 
log(Fx /Fj) =log(Fx) -I- 0.4 J + 6.30, we obtain log{Fx/Fj) = -1.7. 
We note that ROSAT sources con'elating with stellar counterparts extracted 
from either the Simbad or Sloan databas es also show a cutoff in their X-ray 
to infrared flux ratio at or near this Hmit jAgueros et alj2009l) . 



stars were detected in the Einstein GPS, whilst flaring, so such 
sources could represent a not insignificant fraction of the popu- 
lation in our current sample. However, on the basis of the above 
evidence it would seem that even allowing for stellar flaring, the 
saturation limit is applicable to the vast majority of our soft band 
detections. Similar arguments and constraints apply in the case of 
tidally interacting close binaries, such as RSCVn systems, which 
maintain fast rotation throughout their main-sequence lifetime pos- 
sibly into their later evolution (Gudel 2004 ) and which are amongst 
the most luminous coronal emitters, with Lx in some cases ex- 
ceeding 10^^ erg s"^ (e.g., lMakaro"vll2003l) . 

Figure[T2]also shows an additional X-ray constraint at about 3 
kpc, which is the distance at which the line-of-sight column density 
exceeds Nh ~ 1-2 x 10^^ cm~^ (albeit based on a very simplistic 
absorption model as detailed in the figure caption). At this value 
of Nh the soft X-ray absorption is such as to move a source from 
the soft to the hard category (assuming a powerlaw spectrum with 
r = 2.5). The implication is that there is a "spectral" horizon for 
the soft X-ray source sample at about 3 kpc, which corresponds to a 
soft-band X-ray luminosity of somewhat less than 10^^ erg s~ ^ . On 
the basis of the saturation limit for coronal emission, it would seem 
that, at least amongst single stars, only F stars are likely to approach 
the Nh horizon. Conversely, one might conjecture that relatively 
distant F stars are likely to be predominant as one approaches the 
HR=0 boundary. 

Finally, we note that the soft source population could also 
include some accretion-powered stellar systems, in particular 
CVs containing only weakly magnetized white dwarfs. Such 
sources c omprise ~ 90 % of the accreting white dwarf popula- 
tion (e.g., IWamedll995b . have relatively soft thermal components 

s typically 

D^^-^^ erg (e.g.. lVerbuntetalj|l997() . 
non-magnetic CVs might be detected out to the Nh horizon and 
given the low-mass nature of the secondaries be rather faint in the 
near-infrared. However, given the above statistics, this class of ob- 
ject can, at best, represent only a small fraction of the soft source 
population. 

6.2 Nature of the Hard X-ray Source Population 

In the hard X-ray band, if we again take the nominal sensitivity 
limit to be 2 pn ct/ks, then the corresponding unabsorbed flux limit 
is Fx = 2.5 X 10"" erg cm"^ s"^ (2-10 keV). This translates 
to a luminosity of 3 x 10'"' (d/lkpc)^ erg for a source at a 
distance d in kpc. In the 2-10 keV band, our XMM-Newton sur- 
vey is roughly an ord er of magnitude deeper than the ASCA GPS 
dSugizaki et alj|200ll) but shallower by about the same factor with 
re spect to the Chandra d eep pointing in the Galactic Plane reported 
bv lEbisawa et al.l ( l2005h . 

The incidence of NIR counterparts brighter than J=16 within 
the X-ray error circles shows a sharp downward step around HR 
= (Fig. lot, indicative of the fact that the underlying source pop- 
ulations differ markedly between the hard and soft bands. On the 
basis of the 2MASS and UKIDSS cross-correlation analyses, it ap- 
pears that true counterparts may be present, as the brightest star in 
the X-ray error circle, in roughly 25% of cases, down to a limiting 
magnitude of J=20, set against a chance rate which is 2.7 times this 
value. If the cut is set at a brighter limit, at say J < 16, then true 
counterparts may be present in roughly 15% of cases, set against a 
chance rate of somewhat less than twice this value. This rather poor 
"success rate" for an X-ray survey with source locations (or more 
specifically 90% error circle radii) of better than 5" radius is, of 



in their spectra and soft X-ray l uminosities typically in the range 
Lx ~ 1029-32 erg s"^ (e.g.. IVerbunt et alj|l997n . Presumably 
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Figure 17. The limiting distance out to wliicli sources can be detected plot- 
ted as a function of the stellar-spectral class. Here spectral type 0-39 cor- 
responds to dwarf stars from F0V-M9V, with interpolation across spectral 
classifications, as necessary. The thick black lines correspond to the dis- 
tances at which stars reach limiting J magnitudes of 16 and 20. The visual 
absorption. Ay, in the Galactic plane was assumed to increase at the rate 
2 mag/kpc, comparable to that inferred from the UKIDSS GPS survey data 
for {I b)={3l,0) ( see Lucas et af 2008), with the conversion Aj = 0.286 X 
Ay iFiorucci"& M unari 2003). The horizontal grid illustrates the distances 
for a given X-ray luminosity out to which sources may be detected assum- 
ing an effective X-ray survey limit of Fx = 4 X lO"'^^ erg cva~^ 
(0.5-2 keV). The soft X-ray absorption in the Galactic plane was determined 
assuming Nj{ = 1.79 X 10^^ Ay cm~^ and a power-law source spec- 
trum with r = 2.5. The lower (red) curve shows the distance to which 
sources may be detected with the median Lx of nearby st ars as repre- 
sented by the NEXXUS catalogue jSchmitt & Liefkd |2004 . The dotted 
(blue) curve represents the distance out to which stars might be detected 
assuming Lx = lO^'' X Z/toj. The horizontal line labelled "NH limit" 
represents the distance at which, due to absorption, sources switch from a 
soft to a hard characterisation. Absolute Mj and Mj,oi val ues for dwarfs 
stars were taken from Table 5 o f lKraus & HillenbranjlioOTl . 



course, evidence for the relative faintness in the NIR (and optical) 
of the bulk of the true counterparts. For the hard-band sample, the 
average A'^^ inferred from the X-ray spectral analysis is ~ 3x 10^'^ 
cm~^, implying Ay ~ 17 and Aj ~ 5, so it is not surprising that 
the NIR objects linked to the hard sources are typically very much 
fainter than the NIR counterparts of the soft sources, for which the 
impact of line-of-sight absorption is much less severe. 

In broad terms, the various classes of hard X-ray emitting 
source split into systems containing either a high-mass or low-mass 
star (see below). For the stellar-light component, the difference in 
the absolute magnitude will be at least AA/j ~ 5 (e.g., comparing 
AIj for main-sequence stars earlier than B8 with that for dwarfs 
later than KO; Kraus & Hillenbrand 2007), implying a difference in 
log{Fx / Fj) for high-mass and low-mass systems (at a given dis- 
tance emitting at the same Lx) of at least a factor of 2. Even when 
the total systemic emission is considered, any additional contribu- 
tion to the light from accretion dis ks, hotspots, etc ., generally fails 
to bridge this gap (e.g., as in CVs: lAk et"ai]|2007h . It follows that 
there will be a strong bias towards high-mass systems within the 
subset of sources with bright NIR counterparts. Of course, for rel- 
atively high luminosity systems (Lx > 10^^ erg s~^) detectable 
in the hard band out to a few kpc and beyond, the impact of in- 



terstellar absorption at J will also be factor tending to suppress the 
"identification" rat^ll 

Although our knowledge of the hard X-ray source population 
of the Galaxy at intermediate to faint fluxes is very incomplete, due 
to the difficulty in identifying large samples of objects drawn from 
hard X-ray catalogues, we do know at least in qualitative terms 
which classes of source may be present in significant numbers. 
The largest single contribution is likely to come from CVs con- 
taining an accreting magnetic white dwarf (polars or intermediate 
polars) in which there is substantial hard e mission characterised b y 
thermal temperatures of kT 10 keV (e.g., lEzuka & Ishida|[r999l) . 
This class of source which may account for the very high density 
of X- ray sources observed in the Galactic Centre (g.g.. lMuno et al.l 
I2OO3I) and also for a substantial fraction of the hard unresolved X- 
ray emission, known as the Galactic Ridge, observed both near 
the Galactic Centre and also along the inner quadrant of th e 
Galactic Plane (e.g., Sazonov et al. 2006; Revni vtsev et al.ll2006l) . 
LMXB and HMXB systems emitting at relati vely low Lx may 
also be present in the survey. For example. iPfahl et alj ( l2002i) have 
suggested that neutron stars accreting from the winds of main- 
sequence stellar companions might be plentiful in the Galaxy. A 
similar idea was suggested by Willems & Kolb ( I2OO3I) involving 
pre-low-mass X-ray binaries. Howeve r, on the basis of a popula- 
tion synthesis model, iLiu & Li| ( l2006h conclude that neutron star 
LMXB transients in relative quiescence, LMXB with white dwarf 
donors and rotation-powered pulsars may provide an alternative ex- 
planation of the high density of faint X-ray sources seen in the 
Galactic Centre. Our survey which encompasses star-forming re- 
gions both in the Galactic Centre region and in inner spiral arms of 
the Galaxy, may also contain new examples of the highly embed- 
ded Supergiant Fast X-ray Transients (SFXTs) recently discovered 
by Integral (e.g.. .Sguera et al. 2006) and of Be-s tar X-ray binaries, 
including 7 Cas analogs (e.g.,|Motch_et aLll2007l) . 

Relatively hard X-ray emission can also be produced in stellar 
sources without reliance on accretion power. Shocks produced in 
the unstable winds of massive Wolf-Rayet (WR) and 0-supergiant 
stars can generate emission above 2 keV. This hard emission can 
be greatly enhance in WR/OB binaries, syst ems which appear to be 
relatively c ommon in the Galactic Centre jMauerhan etai]|2009l: 
iMauerhari e t al. 2010). Coronal emission in close, tidally inter- 
acting, binaries such as RSCVn-type systems and related classes 
of object (BY Dra-type and Algol-type binaries) m ay also reach 
temperatures extending in to the hard X-ray regime jSwank et al.l 
ll98lUAgrawaletal.lll98ll) . Similarly, in specific locations such as 
the vicinities of dense molecular clouds and star-forming regions, 
protostars and T-Tauri stars may also be prominent as hard X-ray 
sources, by virtue of their enhanced high tempe rature coronal emis- 
sion, particularly duri ng flaring episodes (e.g., IStelzer et al.ll200(]| ; 
iTsuiimoto et al.ll2002l) . Of course, it remains to be demonstrated 
whether any new types of Galactic X-ray source are present within 
current hard X-ray Galactic surveys; for example, isolated neu- 
tron stars and/or blackholes accreting d irectly from the interstellar 
medium jAgol & Kamionkowskill20o3) . 

An illustration of the likely makeup of hard-band selected 
samples is provided by the recent study of iMotch et al. I JioTch 
who carried out an optical follow-up campaign on a sample of 
sources selected from the XMM-Newton Galactic plane survey of 



^ Although the impact of interstellar absorption will be reduced at K, the 
commensurately higher stellar densities mitigates against any clear advan- 
tage in focusing on the longer wavelength band for the current application. 
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iHands et alj ^2004) . Out of 25 hard sources in the follow-up sam- 
ple, there were tentative identifications with four HMXBs and and 
with one likely Wolf-Rayet coUiding-wind binary system. A num- 
ber of the HMXBs had particularly interesting properties, for ex- 
ample the group included a very absorbed likely hyper-luminous 
star in a system bearing similarity to r) Carina, a new/Be star binary 
belonging to the growing class of 7-Cas analogs and a possible 
supergiant X-ray binary of the type recently discovered by Inte- 
gral. The inferred X-ray luminosities for these massive star systems 
was Lx = 10^^"^*erg s-\ There were also probable identifica- 
tions of 3 CVs with relatively faint optical counterparts (including 
a previously undiscovered magnetic system displaying strong X- 
ray variability) and somewhat more tentative identifications of 2 
LMXBs. The detection of several potential active coronal emitters 
above 2 keV also underlines the possibility that active RSCVn bi- 
naries might contribute significantly to hard-band selected samples. 
Unfortunately, the identification of limited samples of object does 
not translate to detail constraints on the space density, scale height 
and luminosity function of the known classes of object and whilst 
our knowledge remains very sketchy, we are unable to give firm 
answers to many important astrophysical questions. 

One important unanswered question is what fraction of the 
hard-band sources might be extragalactic interlopers seen through 
the high absorption column of the Galactic plane? Since we re- 
quire the X-ray sources in our sample to be pointlike (see §3), this 
presumably mitigates against the presence of a significant num- 
ber of clusters of galaxies, leaving active galactic nuclei (AGN) 
as the dominant extragalactic contaminant. From the investigation 
of the column densities in i]4.4l it is clear that AGN would gener- 
ally have HR > 0.8. In an analysis of the 2-10 keV log N - log S 
relation measured in the Galactic plane based on the combination 
of AS CA, Chandra and XMM-Newton observations, .Hands et al., 
( |2004[) demonstrate that at the flux levels encompassed in our cur- 
rent survey, AGN might contribute as much as 50% of the hard 
sou rce population, with a n even higher fraction likely at fainter lev- 
els tebisawa et d]|2005h . Interestingly this rather closely matches 
the fraction of hard sources with HR > 0.8. However, we stress that 
these estimates are very sensitive to the factor assumed for the hard- 
band transmission through the Galactic plane and will vary signif- 
icantly across our survey region; for example. Hands et al. ( 2004) 
assume a transmission factor of 0.68 appropriate to (l,b)=(20,0). 
Further consideration of how the spatial density of the sources in 
the current survey varies as a function of flux and position, and the 
impact of AGN contamination, will be the subject of a subsequent 
paper (Paper III -Warwick et al. in preparation). 



7 CONCLUSIONS 

We have used the 2XMMi source lists pertaining to 116 XMM- 
Newton observations targeted at the inner quadrant of the Galac- 
tic Plane, to construct a sample of serendipitous Galactic X-ray 
sources. The main properties of the 2204 sources which comprise 
the sample are: 

(i) The bulk of the sources have total count rates in the range 
2-100 pn ct ks"^ (0.5-12 keV). In the soft (0.5-2 keV) band, 2 pn ct 
ks~^ corresponds to Fx = 4x10"^^ erg s"^ cm"^ (0.5-2 keV), 
assuming an absorbed power-law spectrum with F = 2.5 and Nh 
= 10^^ cm~'^. The same count rate in the hard (2-12 keV) band 
equates to Fx = 2.5x10"^" erg s"^ cm"^ (2-10 keV), in this 
case assuming F = 1.0 and Nh = 3x 10^^ cm~^. 



(ii) Using a characterisation based on whether the majority of 
the counts were recorded below or above 2 keV, the sample splits 
rather cleanly into 1227 soft sources and 977 hard sources. 

(iii) Both the broad-band hardness ratio (HR) distribution and 
the Band Index (BI) plots reflect a wide spread of underlying 
source spectra. For the soft sources, the X-ray spectra may be rep- 
resented as either power-law continua with F ^ 2.5 or as thermal 
spectra with kT ~ 0.5 keV with Nh ranging from IQ^""^-^ cm~^. 
For the hard sources, a significantly harder continuum form is 
likely, i.e., F ~ 1, with Nh ~ lO'^^^'^'' cm~^. The sources with 
HR>0.8 have column densities commensurate with the total Galac- 
tic line-of-sight value. 

(iv) A high fraction (ij90%) of the soft sources have potential 
NIR (2MASS and/or UKIDSS) counterparts inside their error cir- 
cles, consistent with the dominant soft X-ray source population be- 
ing relatively nearby coronally-active stars. In contrast, the success 
rate in finding likely NIR counterparts to the hard X-ray sample 
is no more than ~ 25% down to J=20, set against a much higher 
chance coincidence rate. The make-up of the hard band population 
in terms likely contributors such as CVs, active binaries, relatively 
quiescent LMXBs/HMXB and other classes of object remains un- 
certain. 

In future papers we will explore the average X-ray spectral 
properties of a subset of relatively bright sources drawn from our 
source sample (Paper II), the log N - log S curves for both the soft 
and hard source samples (Paper III) and the broad-band colours of 
likely counterparts to the X-ray sources (Paper IV). 
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